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Radiation comb generation with extended Josephson junctions 
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We propose the implementation of a Josephson radiation comb generator (JRCG) based on an extended 
Josephson junction subject to a time dependent magnetic field. The junction critical current shows known 
diffraction patterns and determines the position of the critical nodes when it vanishes. When the magnetic flux 
passes through one of such critical nodes, the superconducting phase must undergo a 7T-jump to minimize the 
Josephson energy. Correspondingly a voltage pulse is generated at the extremes of the junction. Under periodic 
driving this allows us to produce a comb-like voltage pulses sequence. In the frequency domain it is possible 
to generate up to hundreds of harmonics of the fundamental driving frequency, thus mimicking the frequency 
comb used in optics and metrology. We discuss several implementations through a rectangular, cylindrical and 
annular junction geometries, allowing us to generate different radiation spectra and to produce an output power 
up to 10 pW at 50 GHz for a driving frequency of 100 MHz. 

PACS numbers: 74.50.+r, 74.81.Fa, 06.20.fb, 04.40.Nr 


I. INTRODUCTION 


The field of optical combs has seen a growing interest in re¬ 
cent years 1,2 . The atomic clocks are extremely stable and have 
sharp resonances; for these reasons, they are used as time and 
frequency reference standards. However, their working range 
is limited to the radio frequency region. This limitation has 
been overtaken only a decade ago. A combination of techni¬ 
cal and conceptual improvements has allowed to extend this 
accuracy to higher frequency up to the optical region. The 
key phenomenon is simple: the atoms are manipulated to emit 
periodic sharp energy pulses which, in the frequency domain, 
correspond to a comb signal with the harmonics of the funda¬ 
mental frequency. Since the generated harmonics show very 
sharp resonances, they can be used as a frequency standard 
in the optical region. The realization of the optical frequency 
comb has paved the way to important applications in optical 
metrology 3 , high precision spectroscopy 4,5 and telecommuni¬ 
cation technologies 1,6 . 

Recently, it has been shown that a similar frequency comb 
can be generated with a dc superconducting quantum interfer¬ 
ence device (SQUID) subject to a time-dependent magnetic 
field 7 . The driving induces 7T— jumps of the superconducting 
phase which are associated to voltage pulses generated at the 
extremes of the device. The voltage pulses sequence in the 
time domain translates, upon Fourier transforming, into a ra¬ 
diation comb in the frequency domain with up to hundreds of 
harmonics of the fundamental driving frequency. 

Here, we show how similar effect can be obtained in an ex¬ 
tended Josephson junction. The underlying physics is similar 
to that discussed in Ref. 7, but the details of the implementa¬ 
tion are different. A setup involving extended junctions opens 
up the possibility for different geometries and, therefore, for 
various power spectra of the emitted radiation. In particu¬ 
lar, since the generated radiation comb structure depends on 
the current-magnetic flux relation of the junction, this latter 
can be properly engineered in order to obtain a desired radia¬ 
tion power spectrum. We discuss the rectangular, cylindrical 
and annular junction designs with their different strengths and 
weaknesses as prototypical examples of extended junctions. 
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FIG. 1. (Color online) (a) Cross section of an extended Josephson 
tunnel junction in the presence of an in-plane magnetic field H gen¬ 
erating a flux piercing the junction in the x direction, t and A 
represent the thickness and London penetration depth of the super¬ 
conductors S , respectively, whereas d is the insulator thickness. Pro¬ 
totypical junctions with rectangular, circular, and annular geometry 
are shown in (b)-(d). L, W, R and r are their geometrical parameters. 

The paper is structured as follows: in Section ?? we intro¬ 
duce the physical mechanism leading to the radiation genera¬ 
tion using extended Josephson junctions. The performance of 
such devices are investigated numerically in Section II, where 
the emitted radiation spectra of junctions with different ge¬ 
ometries are compared. Finally, we gather our conclusions in 
Section III. 

Our system is sketched in Fig. 1(a), and consists of an ex¬ 
tended Josephson tunnel junction composed of two identical 
superconducting electrodes 8-12 . We denote by A and t the 
London penetration depth and the thickness of the supercon¬ 
ductors S , respectively, which satisfy the condition t > A. Fur¬ 
thermore, d is the insulator thickness, whereas tn = 2A + d is 
the magnetic penetration thickness. 

For the sake of clarity we focus on a junction in the short 
limit, i.e., with lateral dimensions much smaller than the 
Josephson penetration depth. In such a case the self-field gen- 




















2 


erated by the Josephson current in the weak-link can be ne¬ 
glected with respect to the externally applied magnetic field 
and no traveling solitons are originated 8,13 . We choose a co¬ 
ordinate system such that the applied magnetic field (77), di¬ 
rected along the x direction, is parallel to a symmetry axis of 
the junction whose electrodes planes lies in the xy plane [see 
Fig. 1(a)]. 

For the sake of presentation, we focus on a rectangular junc¬ 
tion as in Fig. 1(b), keeping in mind that a similar discus¬ 
sion can be extended to junctions with different geometries 
[Figs. l(c)-(d)]. In the limit of short junctions, the approx¬ 
imate behavior of the local phase is <j)(y) = Ky + (p where 
(p 13,14 is the superconducting phase at the center of the junc¬ 
tion, k = 2n<&/ (<FoL) (<f>o — 2 x 10“ 15 Wb is the flux quan¬ 
tum), L is the length of the junction whereas = poHtnL is 
the magnetic flux through the junction and p o is the vacuum 
permeability. By integrating the Josephson current density per 
unit length 13 J c (y) over the junction length we obtain 

, fL/2 

Jj eC ($) = / dy3 c (y) sin= /+ /(Osin (p, (1) 

J-L/2 

where f(t ) = sinc(;r<I>/<I>o), and 7 + is the maximal critical 
current of the junction. The measurable critical current of 
the junction as a function of the magnetic flux is I r c ect (<&) = 
ma X(plj ect (<&). It displays the celebrated Fraunhofer pattern 
which vanishes at the diffraction nodes appearing at d> = n<& o, 
where n is an integer [Fig. 2(a)]. 

The phase jump phenomenon we are interested in can be 
easily described in energetic terms. Assuming that there is no 
bias current, the energy associated to the Josephson current is 

Ej(t) = I I r j ect V{t)dt = — Ejo f(t) cos (p, (2) 

where f(t ) is the same defined above and Ejo = <I>o7+ / (2 n). 
Notice that in the last step of Eq. (2) we have used the second 
Josephson relation <p = (2e/h)V ( t ). Initially f(t = 0) = 1 and 
the minima of the potential energy are found at <jf> = 2k k (with 
k integer). When the magnetic flux reaches the diffraction 
node at = <I>o, Ej vanishes, becoming negative for <f> > <f>o. 
To remain in a minimum energy state, cos (p must change sign, 
which implies that the superconducting phase must undergo a 
7l jump. The original prediction of 7T jumps 8 has also been 
indirectly confirmed via measurements in heat transport ex¬ 
periments performed in temperature-biased Josephson tunnel 
junctions 15,16 . Indeed, as it is explained in Ref. 8, the fact 
that the coherent component of the heat current remains posi¬ 
tive when crossing a diffraction node can be understood only 
if the superconducting phase undergoes a 7T-jump. Possible 
applications of this phenomenon for SQUID devices has been 
extensively discussed in Ref. 7. 

To determine the details of the voltage pulses, such as their 
shape and amplitude, the above discussed energetic picture 
is not sufficient. We rely on the so-called resistively and ca- 
pacitively shunted Josephson junction (RCSJ) model 13,14 in 
which the Josephson junction is modelized as a circuit with a 
capacitor C, a resistor R , and a non-linear (Josephson) induc¬ 
tance Lj arranged in a parallel configuration. We consider a 


sinusoidally-driven magnetic flux with frequency v and am¬ 
plitude £, centered in the first node of the interference pattern 
[Fig. 2(a)], so that 

<b(Q = [1 -£cos(27rw)]. (3) 

As a result, the magnetic flux crosses the nodes of the interfer¬ 
ence pattern at t = (2k + l)/4v, with k integer. Starting from 
the RCSJ model we can write an equation of motion for the 
integrated phase f_ L / 2 dyty (y). Because of the symmetry of 
the problem, this reduces to a RCSJ equation for the phase at 
the center of the junction <p 14 : 

fiC h 

^ ^ <P--4/(0 sin <P = /b- (4) 

We rescale the above equation in terms of adimensional time 
t = 2nvt and, using h/ ( 2e ) = <f>o/27T, we obtain 7 

c ^ + ^-a[f{T)sin ( p-8} = 0, (5) 

where 5 = /#//+, c = 2nRCv and a = 7 + R/(<Fov)- The bias 
current is supposed to be small (5 <C 1) and its effect is to 
impose a preferred direction to the n jumps of the phase. Fur¬ 
thermore, we focus on the limits c <C 1 (overdamped regime) 
and \I+R\ ^ 1, as these two conditions maximize the JRCG 
performance 7 . 


II. NUMERICAL RESULTS 

In the following, all the numerical simulations we discuss 
have been performed for rectangular Nb/AlOx/Nb Joseph¬ 
son junctions subject to an oscillating magnetic flux with fre¬ 
quency v = 100 MHz and amplitude £ = 0.9. As the junction 
parameters we have assumed 17 R = 20 Ohm and 7 + = 0.2 mA, 
while the bias current has been set to 7# = 10 _3 / + . The rela¬ 
tively low frequency we have chosen, besides allowing us to 
achieve an important up-conversion in frequency (see Fig. 4 
and related discussion), also assures that the small capacitance 
of the junction C = 10 fF (corresponding roughly to a junction 
area of 18 ~ 0.2 /im 2 ) has a negligible effect on its dynamics. 

As the critical current crosses the diffraction node at d> = 
T>o, the phase experiences a n jump and a voltage pulse is 
generated across the junction. The shape of the pulse is deter¬ 
mined by the product 7 + R: the larger 7 + R, the sharper the 
voltage pulse. Differently to what happens in the SQUID 
implementation 7 , for a rectangular junction the pulse ampli¬ 
tudes are not the same but show an alternating pattern of lower 
and higher peaks. This is due to the asymmetry in the diffrac¬ 
tion pattern of the critical current near the diffraction node 
[see Fig. 2(a)]. 

This real-time voltage comb could be used in several ways. 
One is as a generator of equally spaced voltage pulses to be 
used in electronics 19 . A second one is as a high precision con¬ 
troller. In fact, because of the Josephson relation, the time 
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FIG. 2. (Color online) Diffraction pattern of the critical current for 
different geometries: (a) rectangular junction, (b) circular junction 
and (c) annular junction with different numbers n of trapped fluxons: 
n= 1 (solid blue line) and n = 2 (dashed green line). The red line 
shows <f>(r)/<Fo oscillating with driving frequency V and amplitude 
£ around different diffraction nodes. 


average of a voltage pulse is actually quantized as a conse¬ 
quence of the 7t jump of the phase: 

2e p *2 

— / dtV(t) = A(p = 7T, (6) 

n Jt x 

where t\ and t 2 are the times in which the jump begins and 
ends, respectively. The phase jumps do not depend on the dy¬ 
namics or on the speed of the node crossing. The only condi¬ 
tion to be satisfied is the crossing of the diffraction node. This 
makes the pulse generation robust against imperfection in the 
dynamics of the junction and the driving. A possible applica¬ 
tion could be the high precise control of a quantum logic gate 
for superconducting based qubits 20,21 . 

The voltage pulse sequence in Fig. 3 has even more interesting 
applications as a radiation generator. In fact, in the frequency 
domain it corresponds to a frequency comb similar to the ones 
used in optics 1 . To test this possible implementation we have 
calculated the power spectrum P vs frequency Q. We first 
compute the Fourier transform of the voltage 

V(Q)= [ T dte in ‘V(t). (7) 

Jo 

The power spectral density (PSD) is then 

PSD(£2) = 2 |V(n)| 2 . ( 8 ) 


Finally, the power P is calculated by integrating the PSD 
around the resonances kv (where v is the monochromatic 
driving frequency) and dividing for a standard load resistance 
Rl of 50 Ohm. This is the power we would measure at a 
given resonance frequency with a bandwidth exceeding the 
linewidth of the resonance. 

To increase the output power, we have considered a linear 
array of N identical junctions connected together via a super¬ 
conducting wire as done for the metrological standard for volt¬ 
age based on the Josephson effect 7,22-24 . The coupling among 
different junctions has been neglected: This condition can be 
realized in practice by a suitable design choice which reduces 
the cross capacitance and the inductance between neighbor 
junctions. In this case the current conservation through any 
z-th junction leads immediately to a set of decoupled RCSJ 
equations of the form (5) 7 . Therefore, the dynamics of the 
junctions are independent and the voltage at the extremes of 
the array is found by summing up the voltages of the sin¬ 
gle junctions. Under the hypothesis that all the junctions in 
the chain generate the same voltage V(t), the total voltage 
drop across the device is simply V tot (t) = NV(t). Accord¬ 
ingly the intrinsic power, that is, the power delivered to an 
ideal load, would scale as N 2 . On the other hand, the extrin¬ 
sic power is less trivial and depends on the detection system 
used 7 . The A 2 -scaling of Josephson junctions performance 
is also a widely studied topic in the context of phase-locked 
Josephson junctions arrays 25-27 . In order to get an estimate of 
the device performance, in the following we have calculated 
the emitted power for arrays of N = 10 3 junctions by dividing 
the total voltage by a standard load resistance Rl = 50 Ohm 
as mentioned above. 

Limitations to this simplified analysis can arise if we must 
take into account the effects of propagation of the emitted 
radiation along the chain 7 . In fact, in our model we have 
considered the device as a lumped element and this assump¬ 
tion breaks down as the length of the chain approaches the 
wavelength of the emitted radiation. As a quantitative es¬ 
timate, the minimum wavelength A m i n (emitted at 50 GHz) 
must satisfy the relation 2A m ^ > L where L is the total length 
of the device 7 . Considering a packing density of the junc¬ 
tions of 5 /im, the above condition is satisfied for N = 10 3 
junctions. Despite being detrimental for the device perfor¬ 
mance, the propagation effects can be taken into account and 
corrected. With a careful design of the device, they could also 
be exploited to amplify the output power at specific work¬ 
ing frequencies. Another limiting factor may be an intrin¬ 
sic property of the device, such as the flux flow through the 
superconductor 28 . This could perturb the external magnetic 
flux that we use to induce the 7T-jumps of the phase. Closely 
related parasitic effects have also been studied recently 29 for 
a similar system based on SQUIDs. Despite being beyond the 
scope of the present work, this remains an interesting issue 
that would require further investigation. 

Figure 4(a) shows the emitted radiation power spec¬ 
trum for a chain of N = 10 3 Nb/AlOx/Nb rectangular 
junctions 14,17,30,31 driven by a 100 MHz oscillating magnetic 
field. As we can see, the device is able to provide a power 
of about 10 p W at 50 GHz (corresponding to the 500-th har- 
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FIG. 3. (Color online) Time dependence of the phase (p (dashed 
line, left axis) and of the voltage V (solid line, right axis) across a 
rectangular Josephson junction. One-directional jumps can be re¬ 
alized by applying a suitable /#. The voltage pulses are generated 
when the interference node is crossed. The parameters chosen for 
the calculations are those typical of a Nb/AlOx/Nb junction 17 with 
R = 20 Ohm, 7+ = 0.2 mA, and we set £ = 0.9. The typical junction 
capacitance is about C = 10 fF and it has been neglected as it does 
not affect the junction dynamics. 

monic of the driving frequency). 

The implementation with extended junctions opens the way 
also for a geometric optimization. Choosing a different junc¬ 
tion geometry affects the critical current of the junction and, 
therefore, the position and the form of the nodes. For the cir¬ 
cular geometry the Josephson current exhibits the known Airy 
diffraction pattern, 



lf rc {<$>) = I- 


J\ (7T<t>/<£>o) . 


(9) 


where J\{y) is the Bessel function of the first kind, <f> = 
liiQHRtn and R is the junction radius. For the annular 
junction 32,33 , the Josephson current takes the form 

irm = h 1 / ' dxxJ n (xw *>/<&(,) sin <p. (10) 

where <f> = Ijl^HRtH, cc = r/R , J n (y) is the nth Bessel func¬ 
tion of integer order, R (r) is the external (internal) radius, 
and n = 0,1,2,... is the number of trapped fluxons in the 
junction barrier. The critical currents for these two junc¬ 
tion geometries are defined as I c c irc (<&) = max^T^ rc (<J>) and 
I“ nn (<£) _ mSLX(p Ij nn , respectively, and they are shown in 
Fig. 2 (b) and (c), respectively. We note that the position of 
the diffraction nodes is different from the rectangular geom¬ 
etry case: The critical currents vanish for non-integer values 
of the ratio d>/<Fo- Moreover, in the annular case, we see that 
the slope of I“ nn (<F) at <f> = 0 differs if the number of fluxons 
n changes. The parameter n can thus be seen as an additional 
degree of freedom which has important effects on the shape of 
the emitted radiation power spectrum [see Figs. 4(c) and (d)]. 
The driving is assumed to have the same periodic behavior, 
oscillating around the diffraction nodes as shown in Fig. 2. 

The power spectrum of the circular junction [Fig. 4(b)] is 
similar to the rectangular junction one. It generates smaller 
output power at high frequency reaching 2 pW at 50 GHz. 
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FIG. 4. (Color online) Power spectrum of the Josephson radiation 
comb generator over a 50 Ohm transmission line. The output power 
is shown in W (left) and dBm (decibel-milliwatt) (right). The cal¬ 
culation is performed for a N = 10 3 chain of Nb/AlOx/Nb extended 
Josephson junctions subject to a V = 100 MHz driving. The cyan 
regions correspond to the insets. The panels refer to different ge¬ 
ometries: (a) rectangular junction, (b) cylindrical junction, (c) annu¬ 
lar junction with one fluxon trapped and (d) annular junction with 
two fluxons trapped. Blue and red points indicate the even and odd 
harmonics, respectively. In Figs, (a) and (c) only the even harmonics 
are present because of the comb-like shape of the voltage pulses. The 
junction parameters were set as those in Fig. 3. 


Particularly relevant is the annular junction case. Here, we 
have an additional controllable parameter: the number n of 
fluxons trapped in the junction. The most interesting situation 
is when there is one fluxon trapped (n = 1) [see Fig. 4(c)]. In 
this case, it is possible to modulate the magnetic flux near the 
vanishing point [see Fig. 2(c)] making the flux driving easier. 
In addition, the diffraction pattern is highly symmetric near 
= 0. This allows one to generate a very precise voltage 
pulse patterns that is eventually reflected in a stronger power 
output at high frequency, as shown in Fig. 4(c) (0.1 nW at 
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50 GHz). 

By varying the number of fluxons, the junction diffraction 
pattern changes [see Fig. 2(c)]: Correspondingly, the dynam¬ 
ics of the junction is different, generating different emitted 
radiation power spectra. Figure 4(d) shows the spectrum gen¬ 
erated by an annular junction chain when two fluxons are 
trapped in each junction. The overall power emitted is smaller 
and the signal is accessible up to ~ 10 GHz. The spectral fea¬ 
tures are very different with respect to the single fluxon ones 
[Fig. 4(c)]. In particular, the lower harmonics (a few multiple 
of v) are now suppressed while the output maximum arises 
around a few GHz. 

The insets in Fig. 4 represent the blow-up of the cyan re¬ 
gions in the main panels: Notice that since the voltage pulse 
signals are almost rectified (Fig. 3), the spectra contain pre¬ 
dominantly the even harmonics, the odd ones being orders of 
magnitude smaller. 

The main sources of error that can limit the device perfor¬ 
mance are the imprecisions in the fabrication process. Small 
differences in the geometry of the junctions, i.e., length L for 
the rectangular junctions and radii R and r for the circular and 
annular junctions, will produce off-sets in the fluxes and de¬ 
lays in the phase jumps. The voltages will still sum up but the 
total voltage pulse shape will be broadened by these effects. In 
the frequency domain, this corresponds to an additional cut¬ 
off at high frequency. Another potential detrimental factor is 
the correction to the dynamics due to the intrinsic junction 
capacitance. However, this effect can be accounted for, mini¬ 
mized or corrected by a proper device design. 

Furthermore, notice that our whole description is done by 
considering the effect of the time-dependent magnetic field 
only, assuming there is no induced electric field on the junc¬ 
tions. It is well known that oscillating magnetic field can gen¬ 
erate in turn electric fields which, especially at high frequency 
(^GHz), may significantly alter the effect we discuss. How¬ 
ever, this problem can be overcome by embedding the junction 
chain inside a suitably designed cavity where the electric (TE) 
and magnetic (TM) modes are spatially separated 34,35 . 

The junction array configurations discussed above are suit¬ 
able for the use as radiation emitters up to 50 GHz. The 
most straightforward way to detect the power in this frequency 
range is to couple the device to a transmission line and to feed 


the signal to a commercial spectrum analyzer. To have access 
to higher frequency we must use different materials (for ex¬ 
ample, YBCO as discussed in Ref. 7) or adopt specific chain 
design. This change must be accompanied with a new de¬ 
tection schemes, for example, by using antennas coupled to 
the device electrodes 7 . Finally, in light of possible implemen¬ 
tations, besides the Nb/AlOx/Nb junctions considered in the 
present work, we signal that other materials could be promis¬ 
ing candidates. For instance Nb/HfTi/Nb junctions 36,37 , be¬ 
ing SNS-like (superconductor-normal metal-superconductor), 
would have the advantage of having almost negligible capaci¬ 
tance, despite having a slightly lower 7 + R product. 

III. CONCLUSIONS 

In summary, we have discussed the possibility to realize a 
Josephson radiation comb generator with extended Josephson 
junctions driven by a time-dependent magnetic field. With a 
linear array of N = 10 3 Nb/AlOx/Nb junctions and a driving 
frequency of 100 MHz, we estimate that substantial power [up 
to ^100 pW] can be generated at 50 GHz (500-th harmon¬ 
ics), opening the way to a number of applications. The de¬ 
vice has room for optimization by modeling the geometry of 
the single junctions, the fabrication materials (see, for exam¬ 
ple, Ref. 7), the driving signal and the array design. The dis¬ 
cussed implementation would have the advantage to be built 
on-chip and integrated in low-temperature superconducting 
microwave electronics. 
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